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Abstract 

One-dimensional repulsive (5-function bose system is studied. By only using the 
Bethe ansatz equation, n-particle partition functions are exactly calculated. From 
this expression for the n-particle partition function, the n-particle cluster integral is 
derived. The results completely agree with those of the thermal Bethe ansatz (TBA). 
This directly proves the validity of the TBA. The theory of partitions and graphs is 
used to simplify the discussion. 
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1 Introduction 



We study statistical mechanics of a one- dimensional gas of Bose particles interacting through 
a repulsive delta function potential. The Hamiltonian for the system with n particles reads 



Throughout the paper, we set h= 1, 2m = 1 and assume the potential is repulsive, c > 0. 

The eigenvalues and the eigenstates are obtained by the Bethe ansatz (BA) method 
1^, |T|, and the quantum inverse scattering method (QISM) ^, |, ^, |, |^ . As a periodic 
boundary condition, the Bethe ansatz (BA) equation is derived. Using the BA equation, 
thermodynamic quantities are calculated by the thermal Bethe ansatz (TBA) |T0| , ^ . In 
TBA, an "interpretation" of the particle density and the state density enables us to define 
the entropy of n-particle system in the thermodynamic limit. 

With the quantum Gelfand-Levitan equation 0, ^ , the field operator is expressed as 
a series of the terms each of which is a product of creation and annihilation operators in 
the scattering data space. The grand partition function is written as a field operator P|. 
Creamer, Thacker and Wilkinson 0] calculated the grand partition function using creation 
and annihilation operators of this system, but the analysis involves a delicate regularization 
that 2tt6{0) is replaced by the volume L. 

To be rigorous, it is desirable to examine these results by a different approach. In this 
paper, we present a method to calculate the thermodynamical quantities only by use of the 
BA equation which is derived exactly as a periodic boundary condition from both QISM 
and BA method. We calculate the n-particle partition function and evaluate the n-particle 
cluster integral. 

The paper is organized as follows. The n-particle partition function and the n-particle 
cluster integral are derived with a method, which we call a direct method, in section 2. This 
method is shown explicitly for n = 3. In section 3, we consider the n-particle case, and 
reformulate the results in graphical expressions. The last section is devoted to concluding 
remarks and discussions. To avoid complexities, the details of a mathematical proof are 
summarized in Appendices A. 



n 




(1.1) 



2 



2 Direct method 



We assume a finite system size L and the periodic boundary condition. It is known that the 



total energy E and the wave numbers kj of the system (lA) are determined by the following 
relations, 



n 



E = J2kl (2.1) 

i=i 

kjL = 2-nmj + ^ A {kj, - kj) , (2.2) 

where rrij are integers or half-integers, and A(/c) is the phase shift of two-particle scattering, 

r Z if AT = odd 
^ W.^i ^fAr ...... (2-3) 



Z + l ifA^ = even, 
rrij < rrij^i, (2.4) 

A (A;) = 2arctan(-), (2.5) 



-TT < A (k) < TT. 

Eq.( p.2|) is called the Bethe ansatz (BA) equation . Main objects to be calculated are the 
ra-particle partition function Z^, 

Z„ = Tre-'^^^ /3 = l/fc^T, (2.6) 

and the n-particle cluster integral 6„; 

" (2.7) 



n>l \n>0 / 



where z = e^^ with the chemical potential /i. By definition, Zq = 1 and we simply have 



h = L-'z, = j 



^e-^'\ (2.8) 

271 



We explain a direct method to evaluate the partition function for 3-particle case. To be 
explicit, the total energy is 



E = kl + kl + kl, (2.9) 



3 



and the BA equation is 



kiL 



2nmi + A (^2 - ki) + A (^3 - ki) 
27rm2 + A (/ci - /^a) + A (/cg - /^a) 
27rm3 + A - k^) + A (A;2 - /ca) 

mi < 7712 < ""^3 G ^- 



(2.10) 
(2.11) 



By use of these relations (|2.9|) ^( 2.11| ), we can calculate the partition function as follows, 

= J2 e-^('=i+'='+fci) 

mi<m2<ms 

6 ^ 2 ^ 3 ^ 

mi,m,2,»n3 mi,m2=ing mi=m2=m3 

= \ Yl j ^^^lt^^2t^rn3e-'^(^'i+^'+A^l)+2-i(rnimi+m2m2+m3m3) 
rhi ,m2 ,^3 

^ ~ / ~ f J 

+ An"e-''(''')+'"'""*" 



- E 



dkidk2dkz 



dm 



dk 



-/3(fcf+fc|+fc|)+iL(fcimi+fe2m2+A:3m3) 



mi ,7712 ,m3 

'fci - ^2 + ic\'^'-'^' fk2-ks + 2c\ '^^-"^ As - A;i + zc^ '^^"^^ 



X 



ki — k2 — ic 



ko - k?. - ic 



ks — ki — ic 



1 





dm' 
dk' 




(K 


-k'^ + icy 
— k'2 — icj 






)+iLk"m" 







7712— 2fhi 



(2.12) 



where \dm/dk\, \dm'/dk'\ and dm" /dk" are the Jacobians to be explained shortly. We 
have written explicitly all the steps of calculations which are common to those for general 



n [[TT], |T^, . In the second equality, we use a symmetry of the BA equation with respect 
to the exchange mj,fcj <-> mii,kii, in the third equality, we apply the Poisson's summation 
formula, and in the last equality, we change variables of integration from m to k. 

The relation between k and m is defined by the BA equation. In ( p.l2| ), k' , m', k" and 
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m" are related as follows, 



k"L 



27im[ + 2A - k[) 

2711712 + ^ (^'1 ~ ^2) 

27rm". 



(2.13) 
(2.14) 



Thus, the Jacobians, \dm/dk\, \dm'/dk'\ and dm"/dk", are given by 
(9m 



(2vr) 



dk 



dm' 



dk' 

dm" 

'dk" 



+ 2L'K {ki - k2) + 2L^K (^2 - k^) + 2L^ir (fcg - ki) 
+3LK {k2 - k^) K (^3 - k^) + 3LK {h - K {k^ - k^) 
+3LK (fci - fcs) K {k2 - fca) , 
+ 3LK {ki - k2) , 

1 



2n 



where 



K{k) = 



dA {k) 
dk 



2c 



k'^ + c^' 



(2.15) 
(2.16) 

(2.17) 



(2.18) 



It is readily shown that all terms except rrij = 0, = and fh" = in ( p.l2| ) exponen- 
tially decay as L gets large. Although the discussion may include these decaying terms, we 
hereafter consider the expressions in the thermodynamic limit. Then, Z3 becomes 

Z, = -J^^^{L' + 6L^K {k, - k2) + 9LK {k^ - k,) K {k, - k,)) e'^'^^+'^i+'^i) 

b J 27T 271 271 

_ i I + SLK {k, - k^)) e-K^'^f +^3^) 



1 

+ 3 



die 

~2^ 



-Le 



-/3(3fc2) 



(2.19) 



It is much easier to show 



2 J 27T 27T ^ ^ 



1 f dk" 
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From (|27^), ( |2.19| ) and ( |2.20|) , the cluster integral 63 is given by 



(2.20) 




1-2-3 = 03 



Fig. 3.1: Vertices of this graph are the elements of Ha, and if x < ?/ then y is drawn "above'' 
X (i.e., with a higher vertical coordinate). It is called the Hasse diagram Jl^ of 



3 r ^ppK{k,-k,)K{h-h)e-^i''^''^'^ 

Z J iTX iTX ZTT 



3 f dk-^ dk<2 



2 J 271 27r 
// 



1 r dk 

+ - 



"i^'"). (2.21) 

Thanks to the effective interaction K{k), 63 consists of only 3 terms. In the next section, we 
show that Zn and 6„ can be calculated in the same way. 

3 Partition function and graphs 
3.1 Partition function 

In order to present a general structure of the partition function, we need to explain some 
terminology in the theory of partition 0]; the partially ordered set (poset, for short). 



Let n(S') denote a set of all partitions of a finite set S. In what follows, [n] means a set 
{1, 2, ■ ■ ■ , n}, and write n„ for n([n]). We define partially order in 11(5') by refinement, that 
is, define x < y & n(S') if every block of x is contained in a block of y. For example, if x G Ilg 
has blocks 137-2-46-58-9 and ?/ G Ilg has blocks 13467-2589, then x < y. Special elements in 
n (S) are Os and l^ such that x > Os and x < Is for all x G n(S'). We write ()„ and 1„ for 



Or„] and lr„,]. Fig.|3.1| shows the partially ordered elements in 113, that is , 1-2-3 < 1-23, 2-13, 



3-12 < 123. In addition, we define Mobius function of n(S') inductively as follows, 
fj.{x,x) = 1 for all X G n(S'), 

fx{x,y) = — iJ,{x,z) for all X < ?/ in n(S'). (3.1) 

x<z<y 
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By use of those definitions, we can show some useful relations. For brevity, M and C 
stand respectively nonnegative integers and complex numbers. (1). Let f,g : 11(5') C, 
then, 



^(a;) = 5^/(1/), forallxen(^) 

y>x 



if and only if 



/ {x)=J2f^ y) 9 {y) , for all a; e n (5) . (3.2) 

y>x 

This is called Mobius inversion formula. (2). Let g, J : 11(5') C and /, /i : N — >• C. If 

EM-);^ = E/W;^ ' (3-3) 

n>0 ■ \n>l ■ / 

h{Ns) = Yl f^{^s,x)g{x), (3.4) 
xmis) 

9{x) = Y.\{J^y)^ (3-5) 

then 

f{n) = J]/i(6„,x)j(x). (3.6) 
We prove this relation in [Appendix A. 



With these two relations, we can derive the partition function Z„ and the cluster integral 

hn in a compact way. First we define /,/;.: N ^ C and g : 11(5) C, 

f{n) = n\Lhn, (3.7) 

h{n) = n\Zn, (3.8) 

h{x) ^ J2 e-^('^?+-+^"), (3.9) 

g{x) = J2 e-'^('='+-+'='0, (3.10) 

where 

xeIl{S), (Ti 3 Sj ^ m[ = nij, x = {(Ti, ■ ■ ■ , (tJ , 5' = {si, ■ ■ ■ , s„} . (3.11) 



Recall that m and k are related by the BA equation. We see that Eq.( p.3|) holds from the 
definition of cluster integral bn, Eq.( p.7|) . From the definition (|3l8|)^( pTT0| ), it is easy to show 
that 



y>x 

/z(6„) = h{n). 
Due to the Mobius inversion formula (|37^) , (|3.12|) is equivalent to 

y>x 

Substituting 0„ for x in this equation, we obtain 



(3.12) 
(3.13) 



(3.14) 



h in) 



X] ^^ (O"' y) 9 (y) 



(3.15) 



This is nothing but the condition 

As we have mentioned in the previous section, we do not include exponentially decaying 
terms. Then, ( |3.10| ) is written as 



dm' 



dk' 



(3.16) 
(3.17) 



The transformation matrix and the Jacobian are given as follows; 



dm' 



J_^S L + E."^.N^"K{k'^„-k'^) if a 
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-N^,K{k'^,-k'^) 



if a ^ a' 



(3.18) 



dk' 



En 

(£U{x)ye^ 



L 



Y.^^J:(u^:^'^''-'^ n \(3.i9) 

<7e?/ /t£V{y)\(Tey /\feGH{i) / 



We explain notations in ( |3.18 ) and ( |3.19| ). The number of elements in a set a is denoted 
by N(j. V (S) denotes a set of all undirected trees with a vertex set S. For instance, all the 
elements in V({(Ti, (J2, cts}) are depicted in Fig.p^ We call a connection of two vertices a 
branch. S (t) denotes a set of all branches contained in a tree t. n{a, t) is the number of 
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Fig. 3.2: A set of undirected trees, V({o"i, o"2, as}), contains three elements 



branches with which the vertex a are connected in the tree t. (Ti{b) and (T2{b) denote two 
end-vertices connected by a branch b. 

A remaining task is a consistency with ( |3.5| ). We define a function J : Il{S) — C by 



J{x) 



HdK 



dm' 



dk' 



(3.20) 



where 



dm' 



dk' 



L 



(2 



Kcrex ) t&V(x) \crex j \f)eH(t) / 



eV(x) \o-gx / \f)eH(t) 

We can change the expression of the Jacobian \dm! jdk'\ from ( p.l9| ) into a sum of forests, 
that is, sets of trees. On the other hand, the right hand side of (|3.21|) is a sum of connected 
forests, that is, trees. Therefore, we put a subscript c as \dm! Idk'\c. With this definition, 
we see that the condition ( p.5|) follows from ( p.l7| ). 

In this way, three conditions (|3.3| )^( ^3D are shown to be satisfied, which indicates that 
Eq.(IJ) holds. We write Eq.Q explicitly, 

dm! 



In fact, (|3.15| ) is equivalent to 



dk' 



dm' 



(3.22) 



dk' 



(3.23) 



The Jacobians in (|3.22| ) and (|3.23| ) are respectively (|3.21|) and (|3.19|) . Explicit form of 
/i(0„,x) can be derived from (^]l]), 

/x(On,x) = \{{-lf^-\N^-l)\. (3.24) 

adx 



Fig. 3.3: Z-toron. 





Fig. 3.4: A set of the trees, V3, contains three elements 

It is easily shown [T^, |T^ that the cluster integrals 6„ agree with those derived from the 
thermal Bethe ansatz (TBA). It is also possible to derive the integral equation in TEA from 

We give two remarks here. First, the Jacobian \dm' /dk'\ is essentially the inner product 
of the Bethe wave functions ||^. Second, the cluster integrals consist of only a finite number 
of terms. 



3.2 Graph representation 

We further develop a graphical representation of the cluster integral and the partition func- 
tion. We draw an I times rolled coil for a Boltzmann weig ht e-^"= (Fig.|3j). We call it 



/-toron following MontroU and Ward [|T5|,|T6[. The tree consists of a toron or torons con- 
nected by branches. The forest consists of the trees. We denote by V„ a set of all the trees 
which satisfy the following two conditions: (1) all the vertices of the tree are made of torons. 



and (2) the sum of rolled number of torons composing the tree is n. Fig. |3.4| shows all the 
elements in V3. W„ is a set of all the forests which satisfy the following two conditions: (1) 
all the vertices of the forest are made of torons, and (2) the sum of rolled number of torons 
composing the forest is n. Simply, is a subset of W^. Fig. |3.5| exhibits all the elements in 
W3. 



In terms of these terminologies, we may rewrite ( p.23[ ) and ( |3.22| ) as 

Sym{f) 



/ew„ 



n 

t in the forest / 



(3.25) 
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Fig. 3.5: A set of the forests, W3. Three graphs in the first hne are not trees because the 
vertices are not connected on partially connected by branch(es). 



Sym (t) 
n\L 



Sit). 



(3.26) 



Here, Sym (t) and Sym (/) indicate symmetrical factors of graphs. In the case t G or 
/ G W^, Sym (t) or Sym (/) means the number of different ways in which a set {1, ■ ■ ■ , n} 
can be distributed to all the vertices of t or / at a time, on the condition that I elements are 
placed in a vertex made of Z-toron. For example 



Sym(^<^\c^ = 12, 
Sym ( = 3, 
Sym {000) = 1. 



(3.27) 
(3.28) 
(3.29) 



From (|3.22| ) and (|3.23|) , we can show that S(t) in ( |3.25|) and ( p.26| ) is given as follows 



X 



dk 



- /C<72(b) 



(3.30) 



,fees(t) 



where denotes the rolled number of toron u, and J2uj (^^ Yluj) denotes a sum (or a 
product) with respect to all the toron u in the tree t. For example. 



S(m) = L X (2 + 1) X ((-l)'l!2° X (-1)"0!1°) 



2ti 2ti 

Substitution of ( p.28| ) and ( |3.31| ) gives the second term in ( p.21| ). 



(3.31) 
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As examples, we list a graphical representation of the cluster integrals 61 ~ 64, 




(3.32) 
(3.33) 



+ QD O + ei) (3.34) 



64 = G +Cf D + c/ u + 

+ & + ® + (3.35) 

4 Concluding remarks and discussions 

In this paper, we have studied a one- dimensional 5-function bose gas, ( [L.lD , and have derived 
directly the partition function and the cluster integral, (|3.23| ) and ( p.22|) , from the Bethe 



ansatz equation, ( |2.2| ). This should be regarded as a proof of the thermal Bethe ansatz 
(TBA). The derivation is simplified by use of the partially ordered set in the theory of 
partition. 

This direct method has some advantages. First, the method is rigorous. In the quantum 
field theoretic method, calculation is done at the infinite volume, therefore an interpretation 
of 27r5(0) as the volume is unavoidable. In the TBA, it is necessary to define the n-particle 
entropy. This procedure is based on an interpretation of the particle and state densities. In 
this sense, both field theoretic method and TBA remained to be proved. On the other hand, 
the direct method is free from such interpretations, and all calculations are traceable step 
by step. 

Second, the method has a wide applicability. It should be remarked that in the direct 
method all the dependencies on the systems come only from the BA equation. In other 
words, the direct method may have the generality at least as well as TBA. 

We conclude that the direct method will be useful and important in calculating the 
thermodynamic quantities in various integrable systems. 
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Appendix A A proof of ( |St6| ) 



We use the same notation as in subsection |3.1| . It can be shown that, on condition that 
h, f : N ^ C, we have 



n>0 



n] 



u 



exp ^ / 



n] 



,n>l 



U 



if and only if 



h{n) = En/(^'^) 
/i(0) = 1. 

This is known as the cumulant expansion formula. Eg. ( |A.2| ) is equivalent to 

\{h{N^) = Y.W^^'^)- 

For H, F : Il{S) C, the following relation holds, 

y<x y^x 



(A.l) 

(A.2) 
(A.3) 

(A.4) 



(A.5) 



This is the Mobius inversion formula, which is a dual form of (|3.2|) . With this formula, 
Eq.( |A.4D becomes 



y<x crGy 



(A.6) 



Substitution of 1„ for x in ( |A.6| ) yields 

fin) = J2fx{y,in)llhiN„). 

We suppose the existence of g : n(S') C which satisfies 



(A.7) 



5; 



xgn(5) 



X] . 



Using this relation, Eq.( |A.6| ) is written as 

fin) - 



xeiin a-ex j/en(cr) 

^fi{6n,x) f^{^,ix)Ylg{y) ■ 

xeu„ 5en(x) ye^ 



(A.9) 
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In the second equality, we have used the following two relations, 

E^(^-')n E ^(y) = E E F{N^)llG{y), (A.IO) 

for F : N ^ C, G : n(5) ^ C, and 

for X G n„. And, we suppose the existence of a map J : n(S') ^ C which satisfies 
Then, a main part of the right hand side of Eq. (|A.9| ) becomes 

E /^(e,yn^(^) = E -"(^'yri e n-^w 

?en(x) ye? Cen(x) ye? CGn(y) zee 

Cen(x) \Aen(C) / zee 

= J(x). (A.13) 



The third equality is due to the following relation. 



J2 f^i^As) = HhNs). (A.14) 

xen(5) 



By use of the relation ( |A.13| ), ( |A.9| ) becomes 



fin)= ^/i(6„,x) J(x), (A.15) 

a;en„ 



which is ( |3.6|) . To summarize, using O, Q) and ( [AI2D , we have proved (| 
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